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DIFFERENTIAL HARNACK INEQUALITIES FOR 
LINEAR PARABOLIC EQUATIONS 

PAUL W.Y. LEE 


Abstract. We prove matrix and scalar differential Harnack in¬ 
equalities for linear parabolic equations on Riemannian and Kahler 
manifolds. 


1. Introduction 

Harnack inequalities are one of the most important estimates in the 
theory of elliptic and parabolic partial differential equations. In the 
linear parabolic case, a version of such estimate was done in [16] in 
connection with the regularity of solutions. For the heat equation, a 
sharp version of the estimate was proved in [T5|. The proof is done by 
proving a gradient estimate, the so-called differential Harnack inequal¬ 
ity. The sharp Harnack inequality is obtained by integrating this gra¬ 
dient estimate along geodesics. Since then, numerous generalisations 
of the differential Harnack inequality to different parabolic equations 
were found (see, for instance, pin El do]). 

The differential Harnack inequalities mentioned above are scalar in¬ 
equalities. A matrix analogue for the heat equation was proved in p. 
A version for the Ricci flow was found in [S] (see P for an application 
and see also PPEldZ] for further developments). 

In this paper, we obtain the following matrix version of differential 
Harnack inequality for general linear parabolic equations on compact 
Riemannian manifolds. 

Theorem 1.1. Let pt be a positive solution of the equation 

pt = Apt + (Vpt, A) + Upt 

on a compact Riemannian manifold. Let A be the operator defined 
by A = VA — (VA)*, where (VA)* is the adjoint of the linear map 
V i-A VA(n). Let W = divX + ^|A|^ — 2Lf. Assume that 

(1) M has non-negative sectional curvature, 

(2) VRc = 0, 

(3) VA = 0, 
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(4) 1^2 ^ y2^ < f^2j 

Then 

vMogp, + i(vx+(v^r)>-^i^^/, 

If X is a gradient vector field, then we recover the following result 
in [TT] , 


Corollary 1.2. Let pt he a positive solution of the equation 


pt = Apt + {Vpt, V/) + Upt 


on a compact Riemannian manifold. Assume that 

(1) M has non-negative sectional curvature, 

(2) VRc = 0, 

(3) v^iy < P/. 


Then 


( log Pt + ^f]> 


k coth(/ct) 


I. 


By setting f = 0, U = 0, and letting A; —)■ 0, we recover the following 
result in [7]. 


Corollary 1.3. Let pt he a positive solution of the equation 

Pt = Apt 

on a compact Riemannian manifold. Assume that 

(1) M has non-negative sectional curvature, 

(2) VRc = 0. 

Then 

V^logpt > 

We also consider the case when M is a Kahler manifold equipped 
with a complex structure J. 


Theorem 1.4. Let pt he a positive solution of the equation 

pt = Apt + {Vpt,X) + Upt 
on a compact Kahler manifold. Assume that 

(1) M has non-negative bisectional curvature, 

(2) J*AJ = A, 

(3) '\/A{X, X) + V A{JX, JX) = 0 for all tangent vector X, 

(4) ^A^ + + J*V^WJ) < k^I. 
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Then 

V' log Pi + log ptJ + + (VX)*) 

+ ^j*{vx + {vxy)j > -kcoth{kt)i. 

The conditions and conclusions are simplified significantly when X 
is a gradient vector field. 

Corollary 1.5. Let pt be a positive solution of the equation 
Pt = Apt + (Vpt, V/) + 17pt 
on a compact Kdhler manifold M. Assume that 

(1) M has non-negative bisectional curvature, 

(2) + J*V^WJ) < k"^!. 

Then 

(^ogPt + + J*X^ (^logpi + J > -kcoth{kt)I. 

By setting f = 0, U = 0, and letting /c —)■ 0, we recover the following 
result in [5]. 

Corollary 1.6. Let pt be a positive solution of the equation 

Pt = ^Pt 

on a compact Kdhler manifold M. Assume that M has non-negative 
bisectional curvature. Then 

VMogpi + J*V‘^ log PtJ > 

L 

Note that the proof of Theorem II .41 is slightly different from that of 
[5]. In particular, the proof here does not make use of any holomorphic 
coordinates. 

There are also scalar versions of the above inequalities with much 
weaker assumptions. 

Theorem 1.7. Let pt be a positive solution of the equation 

Pt = ^Pt + (Vpt, X) + Lf Pt 
on a compact Riemannian manifold. Assume that 

(1) Rc{X,X) > —K\X\'^ for some constant K, 

(2) X^divAl"^ ~ for some positive constants X 

and ki, 

(3) W > —k 2 for some constant /c 2 . 
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Then 


a{t) (^A\ogpt + ^divX^ - b{t) 


Vlogpt + -X 


+ > 


c{t) 


where 


a{t) = sinh^(xt) + {K + A) 


sinh(xt) cosh(xt) 


X 


b(t)^-(K+X) (^'^Hxt)coMxt) _ _ 

c{t) = nsinh(xt) {{K + A) sinh(xt) + xcosh(xt)) 


x = \{K + xy + 


ki + 2(iX + A)/c2 


n 


If X is a gradient vector field, then Al = 0 and we can also set A = 0. 
Therefore, we obtain the following resnlt in m- 

Corollary 1.8. Let pt be a positive solution of the equation 

Pt = Apt + (Vpt, V/) + f/pt 

on a compact Riemannian manifold. Assume that 

(1) Rc > 0 for some constant K, 

(2) AW < ki for some positive constants A and ki. 

Then 


i{t) ^Alogpt + ^A/ 


> 


c{t) 


where 


a{t) = sinh^ | y —t j , 


c{t) = \Jnki sinh / 

On the other hand, if we set X = 0, [/ = 0 in Theorem 11.71 then we 
recover the following result in 


Corollary 1.9. Let pi be a positive solution of the equation 

Pt = Apt 

on a compact Riemannian manifold. Assume that 

(1) Rc{X,X) > —X|Xp for some positive constant K. 
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Then 

a{t)A\ogpt - b{t) iVlogptl^ > 

where 

a{t) = e^^ smh{Kt) — Kt, 

b{t) = — sinh(i^t) cosh(i^t) + Kt, 

c{t) = nKe^^ smh.{Kt). 

By expanding the above inequality in K, we recover the following 
result in jl5j . 

Corollary 1.10. Let pt he a positive solution of the equation 

Pt = ^Pt 

on a compact Riemannian manifold with non-negative Ricci curvature. 
Then 

Tl 

A\ogpt > 
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2. On linearizations of flows 

In this section, we discuss some preliminary results concerning a 
general flow pt of a time-dependent vector held on a Riemannian 
manifold M. 

Let X be a point on the manifold M and let {ni(0), ...,n„(0)} be an 
orthonormal frame in the tangent space T^M. Let Wi{f) be the parallel 
transport of nj(0) along the curve t ha v?t(x) and let K(t) be the skew 
symmetric matrix dehned by 

Kijit) = ^ i{WYt{wi{t)),Wj{t)) - {\/Yt{wj{t)),Wi{t))). 

Then there exists a 1-parameter family of orthogonal matrices 0{t) 
such that 

d{t) = 0{t)K{t). 
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Finally, we define Vi{t) = Oij{t)wj(t) and we call the resulting or- 
thonornial frames {ui(t), parallel adapted frames. Note that 

Vi satishes 


( 2 , 1 ) 


D 

dt 


i=i 
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where £ denotes the covariant derivative, S{t) is the matrix dehned by 

Sij{t) = {VYt{vi{t)),Vj{t))^^ 

and is the skew symmetric part of S{t). 

We remark that similar moving frames adapted to different geometric 
situations were used in the earlier works [mini IS] which, in turn, is 
motivated by a discussion in [T8] . 

For a given matrix A{t), we will denote the symmetric and the skew- 
symmetric parts of A{t) by A^^it) and respectively. Let B{t) 

and R{t) be the matrices dehned by 

= h{Yi + Vy,Fi)(u,(t)),u,(f)\ . 

\ ' 

Lemma 2.1. The matrix satisfies 

s^y{t) + s^y{tfi + s^y{t)s^’^{t) 

- s{t)^\t)s^y{t) + s^^tfi - B^y{t) + R{t) = 0. 


Note that when = V/, we obtain the usual Bochner formula by 
taking the trace of 02.21) . 

Proof. We hrst compute S{t)-. 

Rfit) = (vYfivfit)) + V^Yt{Yt{ipt),vfit)),vfit)'j 

n n 

+ Y^Yt) + (!) (VY,(v.(t)),n(t)) 

k=l k=l 

= (vv,(v.(t)),v,(t)) + (W,(Ki((),r,(v>,)),».,■(()> 

n 

+ (RmK.Mim, »;(<)) + - Sf/(t)Sf*(()) 

k=l 

= (v(Yi + Vr,Yi)(vfit)),vfit))-VYfiVYfivfit))) 

n 

+ (Rm(y;,t,i(i))K„„,(()) + ^(S.f (t)sf/(t) - sf/(t)sf,‘’(i)). 

k=l 
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Therefore, 

S{t) = B{t) - S{tf - R{t) + 

= B{t) - - R{t) - 2S^y{t)S^\t). 

By taking the symmetric part, the result follows. □ 


3. Proof of Theorem 11.11 


In this section, we give the proof of Theorem 11.11 We hrst give a 
proof for the case when the manifold M is compact. 

Let pt be a solution of pt = + (Vpt^X) + Upt and let Yt = 

—2Vlogp4 — A. A computation shows that 


d 

dt 


log Pt = A log Pt + 


V log Pt + -X 


-l\x? + u, 


and 

Yt = VA(-21ogA) - V |-2Vlogpt - A I' 


■\X\^ -2U 


= Vdiv(yi) - V Q \Ytf'^ + V (^div(A) + i|A|2 - 2U 


= Vdiv(y 4 ) - (VYtriYt) + V ( div(A) + -|A|^ - 2U 


Here (VF*)* denotes the adjoint of the linear map v i-A Vht(n). 
Therefore, we have 


(3.1) Yt + Xy^Yt = Vdiv(yi) - A{Yt) + VW, 

where W = divA + i|A|2 - 2U and ^ = VA - (VA)L 
Hence, if B(t) , B(t), F{t) are the matrices dehned by 

Dtj{t) = {VA{vi{t),Vj{t)),Yt ), 

Eij{t) = {X^dW{Yt){vi{t)),Vj{t)) , 

Ftj{t) = {V^W{vt{t)),v,{t)) , 

respectively, then 

Bij{t) = (Xvi{t){Yt + XYtYt),Vj{t)^ 

= Dij{t) + Eij{t) + Fij{t) - {A{Xvi{t)Yt),Vj{t)) 

n 

= + Et,{t) + R^{t) + 2 S,k{t)Sl^{t). 

k=l 
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This gives 


B{t) = D{t) + E{t) + F{t) + + 2{S^\t)f 


and so 


= D{t) + E{t) + F{t) + s^y{t)s^\t) - s^\t)s^y{t) + 2{s^\t)y. 


By combining this with 02.21) . we get 


(3.2) S^y{t) + S^y{tf - - D{t) - E{t) - F{t) + R{t) = 0. 


Next, we consider the term E{t). 

Let us £x a time t and let n be a unit tangent vector which achieves 
the supremum A(t) := sup|^g 7 .^||^|^;^} (V.,,!), n). Assume that v are 
contained in the tangent space T<^j( 3 ,)M. We extend v and Vi{t) to 
vector helds, still denoted by same symbols, dehned in a neighborhood 
of (pt{x) by parallel translation along geodesics. It follows from the 
Ricci identities that the followings hold at 


n 



i=l 


n 




(VFt(ni(t)), Rm(nj(t), n)n) - (L), VRm(ni(t), ni(t), n)n)) j 


n 




{VYt{Rm{vi{t),v)v),Vi{t)) - {VYt{vi{t)),Rm{vi{t),v)v) 
(Yt, VRm{vi{t),Vi{t),v)v)) j 
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n 

= ^ (^{'^^yt{v,Vi{t),v),Vi{t)) - {V^A{v,Vi{t),Vi{t)),v) 

i=l 

- {VYt{Rm{vi{t),v)vi{t)),v) - {VYt{Rm{vi{t),v)v),Vi{t)) 

- {VYt{vi{t)),Rm{vi{t),v)v) - {Yt,VRm{vi{t),Vi{t),v)v))^ 

n 

= (y{'^^Yt{v,v,Vi{t)),Vi{t)) - {V^A{v,Vi{t),Vi{t)),v) 

i=l 

- (VFi(i;),Rm(i;i(t),T)i;i(t)) - (F*, VRm(i;, T)i;i(t))) 

- (VFt(Rm(Ti(t),i;)i;i(t)),i;) - (VFt(Rm(Ti(t), 

- {VYt{vi{t)),Rm{vi{t),v)v) - (Xt,VRm{vi{t),Vi{t),v)v))^. 
Therefore, if Mt = VYt + (VYt)*, then 

n 

A {V,Yt, v) = J2{ V, Vi{t)), v,{t)) 

i=l 

( 3 . 3 ) “ {Mtiv),Rm{vi{t),v)vi{t)) - {Yt,\/Rm{v,Vi{t),v)vi{t))) 

- {Mt{vi{t)),Rm{vi{t),v)v) - {Yt,VRin{vi{t),Vi{t),v)v)) 

- {V^A{v,Vi{t),Vi{t)),v)y 

The parallel adapted frame can be chosen such that {vi{t), ...,Vn(t)} 
form a basis of eigenvectors for the operator (VYt + (VYt)*) with 

eigenvalues Ai,A„, respectively. We can also assume that Ai < ... < 
Xn- Note that A„ = A(f). Once again, we extend Vi{t) to a vector held 
dehned in a neighborhood of ^pt(x) as above. It follows that 

0 > {V^div(Yt){v),v) 

n 

- ^ {V‘^A(v,Vi{t),Vi{t)),v) - (A - Ai) (Rm(T,ni(f))Ti(f),n) 

i=l 

- {Yt,VRm{v,v,Vi{t))vi{t))) + (1), VRm(ni(f), Ti(t), n)n)) j. 

By the Bianchi identity, the above equation becomes 
0 > (V^div(h^)(n), v) + 2VRc(t, v, Yt) — VRc(h), v, v) 

n 

+ ^ (^ - {V^A(u,Vi{t),Vi{t)),v) + (A - Ai) {Rm{v,vi{t))vi{t),v)y 

i=l 

If we assume that the sectional curvature is non-negative and that the 
Ricci tensor and the operator A are parallel, then (V^div(l^)(T), v) < 
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0. Therefore, (13.21) becomes 

X{t) + X{tf - V) - {F{t)V, V) 

= (^s^y{t)v, - s^\tf)v, V) - {F{t)v, V) 

<0, 

where V is the coordinate vector of v relative to the basis 
Finally, if we assume + V^VF < k^I, then 

x{t) + x{ty - e < 0. 

An argument using Gronwall’s inequality shows that 

X{t) < kcoth{kt) 

and the result follows. 

4. Proof of Theorem 11.41 

In this section, we give the proof of Theorem 11.41 
Assume that the dimension of the manifold is given by n = 2N. As 
before, •••,n 2 v(t)} denotes the parallel adapted frame along the 

curve where ^pt is the flow of the vector held F) = —2V log pt — X. 

Let J be the complex structure and let J{t) be the matrix dehned by 
Jiji't) = Assume that A = J*AJ. Since = 

— I {A{vi{t)),Vj{t)), it follows that 

Sf(«) = -i{.4(.,i(()), «,(()> 

.. 2N 

k,l=l 

Since VjL = 0, it follows from fl2.ip that 

j{t) = = 0. 

Therefore, J{t) is independent of time and we can assume that VN+i(t) = 
Jvi{t) for alH = 1,..., TV. It also follows from this and fl3.2l) that 

+ JS^y{t)J^) 
at 

= 

+ + D{t) + JD{t)J^ + E{t) 

+ JE{t)J^ + F{t) + JF{t)J^ - R{t) - JR{t)J^. 


(4.1) 
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Next, recall that the Riemann curvature tensor and the complex 
structure satisfy the following properties on a Kahler manifold: 

(4.2) RmiJX,, JX2)X^ = Rm(Xi, X 2 )X 3 

(4.3) Rm(Xi, X 2 )jrX 3 = :rRm(Xi, ^ 2)^3 

Lemma 4.1. On a Kahler manifold, the following relation on the Rie¬ 
mann curvature tensor holds: 

tr {Y ^ VRm{X, X, Y)Y + VRm{JX, JX, Y)Y) 

= tr{Y^ VRm{Y, Y, X)X + VRm{Y, Y, JX)JX). 

Proof. Let Xi, ...,X 2 n be an orthonormal frame. By (14.21) . (14.dh . and 
the Bianchi identities, 

2N 2N 

VRm(X, X, JXi)JXi + ^ VRuiiJX, JX, JXf)JXi 

i=l i=l 

2N 2N 

= Y, VRm(X, Xi, JX)JX, + ^ VRmiJX, X, Xf)JX, 

2=1 2=1 

2N 

= -Y VRm(Xi, JX, X)JX, 

2=1 

2N 2N 

= Y VRm(X„ JX„ JX)X + ^ VRui{Xi, X, JXi)JX 

2=1 2=1 

27V 27V 

= Y VRm(X„ Xi, X)X + Y VRm(Xi, X^, JX)JX 
2=1 2=1 

□ 

Let n be a tangent vector in We extend n to a vector 

field, still denoted by v, defined in a neighborhood of V5t(x) by parallel 
translation along geodesics. By fl3.3p . fl4.2p . and fl4.3p . 

2N 

A{Vj,Yt,Jv) = {V\dWYt){Jv),Jv)-Y{^"AJv,v^{t),v,{t)),Jv) 

2 = 1 
27V 

+ ^ {J*MJv, Rm(n, Vi(t))vi{t)) + (Yt, VRm(jLn, Jv, Vi{t))vi(t))) 
2=1 

- {J*MJVi{t),Rm{vi{t),v)v) - {Yt,VRm{vi{t),Vi{t),Jv)Jv))y 
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Therefore, by (13.31) and Lemma [4.1[ 

A{{V.Yt,v) + {Vj,Yt,Jv)) = {V\dWYt)iv),v) + {V\dWYt)iJv), Jv) 

2N 

- '^{{'V‘^A{v,Vi{t),Vi{t)),v) + {V‘^A{Jv,Vi{t),Vi{t)),Jv)) 

i=l 

2N 

+ ^ ( {{Mt + J*MtJ)v, Rm(i;, Vi{t))vi{t)) 

i=l 

- ((Mt + J*MtJ)vi(t),Rm(vi(t),v)v} j . 


Let us fix a time t and let Ai < ... < Xn be eigenvalues of M.{(pt{x)) + 
with eigenvectors ni(f),..., i;Ar(f), respectively. Note 
that Jvi{t) is an eigenvector of the same operator with eigenvalue 
Ai. Let A(t) = sup|^g 7 .iv^||„|=i}((VFt(n),n) + (VFi(jLn),jLn)). Then 
A(t) = Xn- We also assume that v{(ftix)) be a unit tangent vector 
which achieves this supremum. It follows that 

0 > {V^{dWYt){v),v) + {V^{dWYt){Jv),Jv) 

2N 

- '^{{V‘^A{v,Vi{t),Vi{t)),v) + {V‘^A{Jv,Vi{t),Vi{t)),Jv)) 

i=l 

N 

+ - ^i) ( (R.m(p, Vi{t))vi{t),v) + (Rm(n, Jvi{t))Jvi{t),v))^ . 

i=l 

Assume that the bisectional curvature is non-negative. Then we 
obtain 

0 > (V2(divFi)(n),n) + {V^{divYt){Jv), Jv) 

2N 

- ^((V^A(v,Vi(t),Vi(t)),v} + (V^A(Jv,Vi(t),Vi(t)),Jv}). 

i=l 

If we assume that '\7A{X,X) -|- VA{JX,JX) = 0 for all tangent 
vectors X, then 


0 > {V\d\wYt){v),v) + {V\dRYt){Jv),Jv) . 
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If we let V be the coordinate vector of v with respect to V 2 N(i)}, 

then fl4.ll) becomes 

A(t) < - V) - V) + ^ {A\v),v) 

+ ^ {A‘^iJv),Jv) + {\/^W{v),v) + {\/^WiJv),Jv) 

< -2a{t)X{t) + 2aitY + ^ {A'^{v),v) + {V‘^W{v),v) + {V^W{Jv),Jv) . 
Therefore, if + VW + J*V^WJ < 2k‘^I, then 
X{t) < —2a{t)X{t) + 2a{tY + 2fc^, 
where a{t) = kcoth.{kt). It follows that A(t) < 2kcoth{kt). 


5. Proof of Theorem 11.71 

In this section, we give a proof of Theorem 11.71 We use the same 
notations here as that of Theorem 11.11 

Let gt = a(t)div(Pt) + — b{t)W. We will apply maximum 

principle on gt. The functional parameters a{t) and b{t) are given by 


a{t) = sinh^(xt) + N 


sinh(xt) cosh(xt) 


X 


b{t) = -N 


sinh(xt) cosh(xt) 
X 


where N = K + X. 

First, we take the trace in fl3.2p . use Rc(X,X) > —K\X\‘^, and 
obtain 


0 = Mt) + - (div.A.y,)^, - + R(t) 

> 5(0 + ^s(t) - h(A + jK. 

n n 4 

- {divAy.)^. - A(divy,)(<^,) - Aiy(^,) - K\Y,\l^, 


where s{t) = div(F^)((pi). 

It follows from this, fl3.ip . a{t) + b(t) > 0, and 


div(Ft) = -2A ( ^ logpi ) = A(divFt) + AVF - AA\Yt\ 
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that gt satisfies 

= aWdiv(yt)((/7t) + a{t)j^dw{Yt){^t) + 

+ h{t) (y, + - b{t)W{^t) - b{t) {VW,Y,)^^ 

= a(t)div(rt)((^t) + {a{t) + b{t))j^dW{Yt){ipt) + ^\Yt\l^ 

- b{t)dW{Y,){^,) - b{t) {AiY,),Y,)^^ - b{t)W{^,) 

< (aW - div(y,)(^,) + (^ic(t)^ - \\A\l^ (ait) + bit)) 

+ ((div^, Yt)^^ + A(divyt)((/?t) + AWi^pt) + /tr|yt|^J(a(t) + bit)) 

+ ^-^\Yt\l - bit)dWiY,)i^,) - bit) {AiY,),Y,)^^ - bit)Wig^,) 

= (a(t) - div(y,)M + ^giv^t) + ^-^\Yt\l - bit)Wip>,) 

+ ^(div^,+ K\Yt\‘^^^ + -c(t)^ - + AWiipt^ (a(t) + bit)) 

< (^a(t) - div(yt)((/7t) + Agi<^t) + ^\Yt\l, - bit)Wiipt) 

+ + i^ + {Ytllt + + AWiift)^ (a(t) + bit)) 

Let c(t) be the function defined by a(t) — ^db(«(b+fe(b) ^ g_ gj^ce 
bit) + 2(JL + A)(a(t) + bit)) = 0, we also have 

J^gti^t) ^ ^gi^t) - hit)Wi^t) 

^ ^ ~ ^ (aW + K^))- 

Since ^\dWA\‘^ — ^\A\‘^ + AW < ki, b < 0, and W > —k 2 , it follows 
that 

gt + (Vg't, Zt) < —I" (®(^) + ^(^)) + ^gt ~ bit)k2 

+ fci + 2(JL + A)A:2 J (Q-(t) + bit)) + Agt- 


n 
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By maximum principle, we obtain 

9t ^ J ki + 2{K + X)k^ (®('7‘) + b{T))dT 

and the result follows. 
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